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Abstract 

Unsteady  interactions  of  strong  concentrated 
vortices,  distributed  gusts,  and  sharp-edged  gusts 
with  stationary  airfoils  have  been  analyzed  in 
two-dimensional  transonic  flow.  A  simple  and 
efficient  method  for  introducing  such  vortical 
disturbances  has  been  implemented  in  numerical 
codes  that  range  from  inviscid  transonic  small- 
disturbance  to  thin-layer  Navler  Stokes.  The 
numerical  results  demonstrate  the  large  distor¬ 
tions  in  the  overall  flow  field  and  in  the  surface 
air  loads  that  are  produced  by  various  vortical 
interactions.  The  results  of  the  different  codes 
are  in  excellent  qualitative  agreement,  but,  as 
might  be  expected,  the  transonic  small-disturbance 
calculations  are  deficient  in  the  important  region 
near  the  leading  edge. 


Introduction 


Vortical  disturbances  can  have  important 
effects  on  the  airloads  and  the  aeroacoustics  of 
a  variety  of  aerodynamic  devices,  such  as  maneu¬ 
vering  aircraft  and  missiles,  helicopter  rotor 
blades,  and  turbomachinery.  Although  the  numeri¬ 
cal  analysis  of  strong  vortical  flows  has  received 
less  attention  than  the  treatment  of  shock  waves 
over  the  past  decade,  compressible  flow  fields 
with  embedded  regions  of  distributed  or  concen¬ 
trated  vorticity  appear  certain  to  receive 
increasing  attention  in  the  future. 

Several  methods  have  already  begun  to  appear 
in  the  literature  for  treating  concentrated  vor- 
ticit.y  in  finite-difference  computations.  Direct 
approaches,  such  as  Rizzi  and  Erickson,1  Fuji!  and 

p  i  a 

Kutler,  Krause, J  and  Srtnivasan  and  Steger  may 
be  classified  as  vortex-capturing,  by  analogy  with 
the  familiar  shock-capturing  methods  that  are  used 
in  many  computational  fluid  dynamics  (CFD)  codes. 
As  in  shock  capturing,  the  details  of  the  actual 
phenomenon  are  spread  over  several  grid  points. 
Consequently,  the  solutions  in  these  high-gradient 
regions  are  artificially  grid-dependent  and  are 
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susceptible  to  the  effects  of  numerical  dissipa¬ 
tion.  In  the  case  of  concentrated-tip  vortices, 
for  example,  numerical  dissipation  may  destroy  the 
simulation  of  the  core  structure  faster  than  would 

ii 

physical  dissipation,  unless  many  grid  points  are 
clustered  in  the  vortical  region. 

Two  alternate  methods  for  modeling  vortex 
flows  are  illustrated  schematically  in  Fig.  1.  By 
analogy  with  shock-fitting  in  transonic  flows, 
these  methods  may  be  thought  of  as  vortex-fitting. 
The  upper  part  of  the  figure  portrays  the  method 
introduced  by  Caradonna  et  al.,6  and  used  by 
George  and  Chang,6  for  example.  This  technique  of 
one  or  more  branch  cuts  permits  concentrated  vor¬ 
tices  to  be  introduced  into  potential-flow 
formulations. 


The  prescribed-vortex  method,  which  is  the 
subject  of  this  paper,  is  shown  in  the  lower  half 
of  Fig.  1,  The  essential  ideas  are  due  to 
Steinhoff7  and  Buning  and  Steger,8  and  the  tech¬ 
nique  has  been  used  by  the  authors^  11  for  poten¬ 
tial,  Euler,  and  thin-layer  Navier-Stokes  analyses 
of  two-dimensional,  unsteady  airfoil-vortex  inter¬ 
actions.  In  these  applications,  the  structure  of 
the  vortex  is  prescribed,  but  its  path  in  space 
can  be  allowed  to  develop  a3  part  of  the  solution. 
The  prescribed-disturbance  method  is  also  well 
suited  to  and  has  been  used  for  more  general  vor¬ 
tical  flows,  such  as  gusts  with  distributed 
1  ? 

vorticity. 


Many  of  the  interesting  practical  applica¬ 
tions,  such  a3  blade-vortex  interactions  of  heli¬ 
copter  rotors  (e.g.,  Fig.  2),  and  gust  response  of 
aircraft  wings,  are  three-dimensional  unsteady 
problems  that  will  require  enormous  computational 
resources  for  transonic  cases.  Therefore,  it  is 
useful  to  study  the  numerical  methods,  as  well  as 
the  general  nature  of  the  phenomena,  in  two  dimen¬ 
sions,  as  indicated  in  the  figure.  It  is  also 
important  to  establish  the  minimum  level  of  com¬ 
plexity  in  the  governing  equations  that  will  suf¬ 
fice.  In  the  following  sections,  we  review  our 
transonic  small-disturbance,  Euler,  and  thin-layer 
Navier-Stokes  analyses,  especially  the  two- 
dimensional  Implementation  of  the  prescribed- 
disturbance  approach.  Several  representative 
examples  are  given  that  illustrate  the  advantages, 
capabilities,  and  limitations  of  these  different 
formulations,  in  anticipation  of  future  extensions 
to  real,  three-dimensional  problems. 


WAKE  SHEET 


BRANCH  CUT 
TO  DEFINE 
VORTEX.  [0]  =  r 


(a)  POTENTIAL  -  FLOW  VORTEX  FITTING 


q  =  U^  +  qv  +  Aqw 


d% 


jr 

4  n 


cos  ds 


ry  CALCULATED  (FORCE  FREE  PATH) 


(b)  PRESCRIBED  -  DISTURBANCE  VORTEX  FITTING 

Fig.  1  Methods  for  modeling  vortex  flows. 
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Fig.  2  Two-  and  three-dimensional  unsteady  vortex  Interactions. 


Numerical  Formulations 
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The  essence  of  the  prescrl bed-dl 3turbance  or 

Q 

perturbation  method  is  simply  that  each  of  the 
dependent  flow  variables  Is  spilt  Into  a  pre¬ 
scribed  part,  which  defines  (or  Is  consistent 
with)  the  Imposed  vortical  disturbance,  and  a 
remaining  part,  which  Is  obtained  from  the  solu¬ 


tion  of  the  governing  equation  set.  Even  though 
the  governing  equatlon(s)  may  be  nonlinear  and 
Independent  solutions  are  therefore  not  super- 
posable,  we  may  still  decompose  the  dependent 
variables  as  follows: 
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(q  -  qJ 


and  qy  represents  the  specified  vortical  distur¬ 
bance  superimposed  on  a  uniform  free  stream. 

The  following  subsections  summarize  the 
development  of  the  governing  equations  for 
(q  -  qv),  and  present  the  methods  for  solving  them 
numerically.  Although  the  overall  concept  is 
essentially  the  same  for  the  Euler  or  Navier- 
Stokes  formulations  as  for  the  transonic  small- 
disturbance  (potential)  formulation,  the  actual 
derivations  and  solution  techniques  are  rather 
different;  therefore,  they  are  treated  separately. 
The  complete  details  of  each  are  given  in  Refs.  10 
and  11  for  the  former,  and  in  Refs.  9  and  12  for 
the  latter. 

Euler  and  Thin-Layer  Navlei — Stokes  Equations 
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Here  U  and  V  are  the  contravariant  velocities 
along  the  £  and  n  directions,  i.e., 

U  -  Et  ♦  ?xu  ♦  5yv  (5a) 

V  -  nt  ♦  nxd  ♦  nyv  (5b) 

and  J  is  the  transformation  Jacobian.13 

In  Eq.  (2)  qv  represents  any  prescribed 
nonuniformtty  superimposed  on  an  otherwise  uniform 
free  stream.  In  the  present  Investigation,  a 
concentrated  vortex  is  considered,  either  fixed  in 
space  or  moving  with  the  flow.  The  velocity  field 
of  the  vortex  is  specified  as 


In  this  part  of  the  investigation,  we  use  the 
unsteady,  two-dimensional  equations  in  strong 
conservation-law  form , 13  including  the  Baldwin- 
Lomax1*4  algebraic  eddy-viscosity  model  for  thin- 
layer  turbulent  flows.  In  generalized  coordinates 
(•;,n,x)  and  in  the  perturbation  form  of  Buning  and 
Steger,8  the  equations  are 

3 t ( q  -  qv)  ♦  3j.(E  -  Ey)  ♦  3n(F  -  Fy)  -  Re'^S 


where  the  dependent  flow-field  vector  is 


and  the  flux  vectors  E,  F,  and  S  are  given  by 
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where^  is  the  dimensionless  vortex  strength 
and  r(x,t)  is  the  instantaneous  radial  distance 
from  the  vortex  center.  That  is,  the  vortex- 
induced  field  qv  is  irrotational  outside  a  vis¬ 
cous  core  of  radius  a.  Thus  the  vector  compo¬ 
nents  p,  (pu),  (pv),  and  e  are  determined  from 
the  inviscid  (Euler)  components  of  Eq.  (2),  as 
discussed  in  Refs.  10  and  11. 

The  usual  tangency  and  no-slip  boundary  con¬ 
ditions  on  the  airfoil  are  applied  at  each  instant 
to  the  total  flow-field  vector  q.  Surface- 
conforming  C-type  grids  are  used,  generated  by  an 
algebraic  method,18  with  clustering  near  the  lead¬ 
ing  and  trailing  edges  and  in  the  vicinity  of  the 
upper  and  lower  surface  shock  waves.  Typical 
grids  for  viscous  calculations  consist  of  161  to 
1 8 1  points  in  the  surface,  £,  direction,  and 
52  points  in  the  normal,  n,  direction,  extending 
6  to  10  chords  in  all  directions.  Forty-five  grid 
points  in  the  normal  direction  are  generally  used 
for  the  inviscid  cases. 

An  implicit,  spatially  factored  algorithm 
with  Euler-implicit  time-differencing1 3  is  used  to 
solve  Eq.  (2).  The  scheme  is  first-order  accurate 
in  time  and  second-order  accurate  in  3pace.  Sec¬ 
ond-order  implicit  smoothing  and  fourth-order 
explicit  dissipation  terms  are  added  to  improve 
the  nonlinear  stability.  Even  so,  the  nondimen- 
slonal  time  steps  generally  have  to  be  restricted 
to  values  on  the  order  of  0.005  during  that  por¬ 
tion  of  the  Interaction  when  the  vortex  is  close 
to  the  airfoil.  The  CPU  solution  time  on  the  Ames 
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Cray  X-MP  computer  la  approximately  2  sec/time 
step.  The  tnvlacid  cases  considered  In  this  paper 
required  approximately  10-15  min  each,  and  the 
viscous  calculations  about  25*30  min  CPU  time, 
after  the  Initial,  steady  solution  was  obtained 
without  the  vortex. 


however,  the  continuity  and  momentum  equations  can 
be  combined  to  give12 
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where 


The  effectiveness  of  the  present  prescribed- 
disturbance  scheme  Is  illustrated  in  Fig.  3  for 
one  of  the  examples  to  be  discussed  in  more  detail 
In  the  section  on  Unsteady  Vortex  Interactions. 

The  figure  shows  the  variation  In  lift  on  the 
airfoil  as  a  strong,  concentrated  vortex  passes 
beneath  it,  with  and  without  the  special  treatment 
outlined  above.  The  "nonperturbation"  scheme 
essentially  represents  vortex  capturing;  that  is, 
the  vortex  structure  Is  specified  only  as  an  ini¬ 
tial  condition,  when  the  center  of  the  vortex  is 
5  chords  up3tream  of  the  leading  edge.  Subse¬ 
quently,  this  Initial  disturbance  convects  through 
the  computational  domain,  and  the  numerical  dissi¬ 
pation  associated  with  the  finite  grid  spacing 
progressively  weakens  the  grad  ents  and  reduces 
the  effective  vortex  strength  substantially.  This 
numerical  error  is  clearly  grid-dependent;  how¬ 
ever,  it  is  completely  absent  in  the  prescrlbed- 
iisturbance  solutions,  which  are  essentially  inde¬ 
pendent  of  the  grid. 


181  x  45  GRID 
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Mg.  3  Effectiveness  of  the  prescribed- 
disturbance  method  for  Euler  calculations. 


Transonic  Small-Disturbance  Potential 
Equat ion 
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dr  is  an  incremental  position  vector  which  i3 
integrated  over  a  suitable  path  from  some  upstream 
reference  point,  and  the  symbol  q  now  represents 
the  velocity  field  (u,v). 

We  now  go  a  step  beyond  Eq.  (1)  and  decompose 
this  total  velocity  field  into  the  following  com¬ 
ponents;  1)  the  uniform  free  stream,  U^;  2)  the 
prescribed  disturbance  field,  q  ;  and  3)  an  lrro- 
tatlonal  perturbation  velocity  ifield  V$  produced 
by  the  airfoil  In  the  presence  of  the  concentrated 
vortex  or  distributed  gust.  That  is, 

q  *  U^  *  qv  ♦  V*  (9) 

Then,  for  example,  for  a  distributed  vortical 
(rotational)  gust 

Q  =J(aT  *  q  x  V  x  qvj  •  dr  ♦  |f-  dOa) 


whereas  for  an  irrotational  moving  vortex 
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at 


(10b) 


where  $y  »  2mr  tan  '[(y  -  y^/tx  -  xt)],  and 
(xj.y,)  is  the  instantaneous  position  of  the 
center  of  the  concentrated  vortex. 


It  should  be  emphasized  again  that  the  decom¬ 
position  of  the  velocity  in  Eqs.  (1)  or  (9)  does 
not  imply  linearity.  Both  the  boundary  conditions 
and  the  governing  equation  for  the  airfoil  distur¬ 
bance  potential,  7$,  are  altered  by  the  introduc¬ 
tion  of  qy,  and  independent  solutions  are  not 
superposable  for  transonic  flows. 


In  this  case,  the  approximations  are  a  thin 
lirfoil  in  an  inviscid,  l3entropic  fluid,  with 
either  (1)  a  concentrated  potential  (irrotational) 
vortex,  or  (2)  small  distributed  vortical  distur¬ 
bances,  superimposed  on  an  otherwise  uniform, 
nearly  sonic  free  stream.  As  in  the  above  formu¬ 
lation,  these  prescribed  vortical  disturbances  are 
assumed  to  not  be  altered  by  the  Interaction  with 
the  airfoil. 


The  governing  equation  for  ^  is  obtained  by 
substituting  first  the  quantity  q  *  U^  *  q  ,  and 
then  q  given  by  Eq.  (9),  into  Eq.  (7).  Next, 
the  former  resulting  equation  is  subtracted  from 
the  latter.  Then  the  quantity 

J(aT  -  do  x  7  x  %)  •  d; 


The  combination  of  unsteady  flow  and  rota¬ 
tional  disturbances  in  the  free  stream  requires 
some  special  attention  if  potential-flow  concepts 
are  to  be  retained.  Under  the  previously  men¬ 
tioned  lsentropic,  small-disturbance  assumptions. 


which  is  common  to  both,  drops  out.  Even  if 
V  x  qy  •  0,  the  quantity  (7$  x  7  x  qy)  -dr 
vanishes  if  dr  is  aligned  with  7$.  Finally, 
the  local  speed  of  sound,  a,  is  evaluated  using 
the  lsentropic  relations  and  the  usual  tran 
sonic  small-disturbance  scaling  laws  and 
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approximations,1®"17  coupled  with  the  observation 


that  the  prescribed  disturbances  considered  In 
this  paper  produce  only  second-order  variations  In 
density,  pressure,  and  temperature.  This  gives 
the  following  modified  form  of  the  unsteady 
transonic  3mall-disturbance  equation  (In  strong 
conservation  form) 


M^tt  *  2M^xt  ‘  b  CCb2  *  C2*x)(*x  *  Uv> 
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Linear  theory  can  be  simulated  by  setting  C2  *  0. 


The  boundary  condition  on  the  airfoil  is 
obtained  by  simply  substituting  Eq.  (9)  into  the 
small-disturbance  approximation  to  flow  tangency, 
gi ving 


V  3y )  Ut  «  3x/ 


T.  -  v 
b  v 


(12) 


As  noted  above,  neither  small -ampl itude  gusts 
nor  moving  concentrated  potential  vortices  produce 
first-order  pressure  fluctuations  in  the  free 
stream.  Therefore,  consistent  with  our  basic 
assumption  that  the  prescribed  vortical  distur¬ 
bance  remains  undistorted  as  it  passes  by  the 
airfoil,  we  retain  the  usual  form  for  the  small- 
disturbance  pressure  coefficient, 


-2(4 


.) 


(13) 


Numerical  experimentation  revealed  that  the 
most  serious  deficiencies  of  the  small-disturbance 
formulation  are  1)  a  limitation  on  the  strength  of 
the  vortex  3nd/or  its  closest  proximity  to  the 
airfoil,  and  2)  the  leading-edge  behavior  for 
strong  interactions  with  very  thin  airfoils.  The 
latter  can  be  partially  overcome  with  an  ad  hoc 
leading-edge  correction  to  Eq .  (12);  viz. 


(in  . 
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as  discussed  in  Ref.  10.  This  treatment  has  been 
found  to  be  very  satisfactory  for  thin  airfoils, 
such  as  the  NACA  64A006  considered  in  the  section 
on  Unsteady  Vortex  Interactions.  However,  compar¬ 
isons  with  Euler  and  full-potential  solutions  show 
that  It  over-corrects  for  moderately  thick  air¬ 
foils.  such  as  the  NACA  0012,  and  that  It  Is  com¬ 
pletely  unsatisfactory  for  airfoils  with  blunter 
leading  edges. 


Equation  (’1)  is  solved  using  the  alternat- 
ing-dlrectlon-lmpl tcit  numerical  algorithm  of  the 


well-exercised  transonic  code  TRAN2 
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of  the  high-frequency  terms  Included 
sions  added  for  variable  time  steps. 


wl th  all 
and  provl- 
A  rectangu¬ 


lar  cartesian  grid  (113  or  186  x  97)  la  used,  with 
clustering  near  tne  leading  and  trailing  edges. 

The  CPU  solution  times  on  the  Ames  Cray  X-MP  com¬ 
puter  range  from  approximately  1  to  A  min. 

Comparison  of  Initial  Solutions 

A  preview  of  the  similarities  and  differences 
in  representati ve  initial  solutions,  representing 
the  steady-state  solutions  of  airfoils  in  uniform 
free  streams,  Fig.  4,  is  helpful  before  consider¬ 
ing  the  vortical  interactions  in  the  section  on 
Unsteady  Vortex  Interactions.  In  this  and  subse¬ 
quent  figures,  the  Reynolds  number  is  6  x  10®  for 
the  thin-layer  Navier-Stokes  solution,  with  a 
turbulent  boundary  layer  assumed  for  the  entire 
airfoil.  The  leading-edge  correction  in  Eq.  (12b) 
was  used  for  the  small-disturbance  calculations  on 
the  NACA  64A006  airfoil,  but  not  for  the 
NACA  0012.  The  small-disturbance  deficiency  in 
the  leading-edge  region  is  more  evident  in  the 
results  at  Mb  «  0.3  than  at  higher  Mach  number. 

For  the  NACA  0012  airfoil  at  -  0.8,  the 
shock  wave  is  fairly  strong,  but  viscous  effects 
are  relatively  weak.  Therefore,  the  Navier-Stokes 
and  Euler  solutions  are  approximately  the  same, 
but  the  details  near  the  shock  wave  are  slightly 
different  in  the  small-disturbance  solution.  On 
the  other  hand,  the  viscous  effects  are  stronger 
on  the  64A006  airfoil,  but  the  Euler  and  small- 

disturbance  are  almost  identical.  The  experimen- 

1 8 

tal  results  are  from  Zwaan. 


Example  Results--Unsteady  Vortex  Interactions 

For  the  sake  of  brevity,  we  consider  only 
concentrated  vortices  convecting  with  the  flow 
past  symmetrical  airfoils  at  zero  angle  of  attack. 
Therefore,  all  the  differences  between  the  flow 
fields  above  and  below  the  airfoils  are  due  solely 
to  the  vortex  interaction.  Examples  of  vortices 
fixed  in  space  and  of  airfoils  at  Incidence  are 
given  in  Refs.  9-1 1  . 

The  small-disturbance  results  in  Fig.  5  illu¬ 
strate  many  of  the  essential  features  of  transonic 
airfoi 1-vortex  interactions.  For  this  case  of  a 
vortex  with  clockwise  circulation  passing  beneatn 
the  airfoil,  the  vortex-induced  "downwash"  effect 
first  appears  as  an  effective  negative  angle  of 
attack.  Figs.  5a  and  5d.  Unlike  the  linear 
results,  however,  the  transonic  effects  are  much 
larger  on  the  lower  surface  than  they  are  on  the 
upper  surface.  Also,  the  actual  shape  of  the 
pressure  distributions,  the  shock-wave  strengths 
and  positions,  and  the  time-history  of  the  airfoil 
pressures  are  quite  different,  even  qual i tat  1 vely , 
from  the  results  for  airfoil  oscillations,  or  from 
the  results  for  the  sinusoidal  gusts  considered  in 
the  section  on  Airfoil-Gust  Interactions. 

As  the  vortex  passes  beneath  the  airfoil, 
Figs.  5b  and  5e,  the  u-component  of  the  vortex- 
induced  velocity  field  also  becomes  important  in 
the  nonlinear  case.  After  the  vortex  passes 


THIN  LAYER 
NAVIER  STOKES 


EULER 


SMALL  DISTURBANCE 


NACA  0012  AIRFOIL,  M^  =  0.8,  n  =  0J 


Fig.  il  Initial  solutions  without  a  vortex. 


b^mnd  the  trailing  edge.  Figs.  5c  and  5f ,  its  and  the  return  of  the  flow  around  the  airfoil  to 

induced  field  becomes  and  effective  "upwash"  its  orginal  state  is  extremely  slow.  The  moving- 

distribution  of  vertical  velocity,  and  its  vortex  interaction  is  substantially  different  from 

horizontal  component  becomes  relatively  unlm-  the  quasi-steady  response  that  would  be  obtained 

portant  again.  Throughout  the  interaction,  non-  by  freezing  the  vortex  at  different  x-positions. 

linear  unsteady  lag  effects  are  very  important, 


Fig.  5  Transonic  airfoil-vortex  interaction  on  the  NACA  0012  airfoil  at 

M  =  0.8,  a  -  0,  r  -  0.20,  y  -  -0.26. 

“  v  ’  Jv. 


Figures  6  and  7  compare  the  transonic  small- 
disturbance,  Euler,  and  thin-layer  Navier-Stokes 
results  for  the  NACA  6 U A00 6  airfoil  at  M  -  0.85, 

0D 

with  the  same  vortex  strength  and  initial  position 
as  before.  This  airroil,  which  has  been  used  in 
numerous  numerical  and  experimental  studies,  is 
not  only  thinner  than  the  NACA  0012  section,  but 
has  a  significantly  smaller  leading-edge  radius  as 
well.  Therefore,  it  might  be  expected  to  be  more 
sensitive  to  the  vortex-induced  downwash  than  the 
NACA  0012  profile. 

As  in  the  steady  results  shown  in  Fig.  4,  the 
3hock  wave  in  this  case  is  smeared  by  the  pre¬ 
dicted  vlscous-invlscid  interaction,  according  to 
the  Navier-Stoke3  calculations.  However,  no 
boundary-layer  separation  is  Induced.  Otherwise, 
the  predictions  of  all  three  methods  are  similar, 
and  the  sequence  of  events  is  essentially  the  same 
as  those  discussed  in  connection  with  Fig.  5. 
However,  larger  fluctuations  appear  in  the  pres¬ 
sure  distributions  near  the  leading  edge.  Without 
the  leading-edge  correction  of  Eq.  (12b)  this 


effect  was  over-predicted  by  the  small-disturbance 
calculations. 

For  the  results  shown  in  Figs.  6  and  7,  the 
vortex  is  constrained  to  convect  at  along  a 

straight  line,  yv  -  -0.26.  Figure  8  shows  the 
effect  of  allowing  the  vortex  to  move  along  a 
force-free  path;  that  is,  along  a  streamline.  The 
Euler  calculations  show  some  differences  caused  by 
this  effect  as  the  vortex  passes  near  and  beyond 
the  trailing  edge,  but  this  behavior  is  not 
evident  in  the  small-disturbance  results. 

Instabilities  usually  develop  in  the  small- 
disturbance  calculations  when  cases  with  sig¬ 
nificantly  stronger  vortex  interaction  are 
attempted.  The  Euler  and  Navier-Stokes  calcu¬ 
lations  require  a  reduction  in  the  nondimenslonal 
time-step  and  require  special  attention  to  the 
dissipation  and  smoothing  of  parameters,  however, 
solutions  have  been  obtained10,11  for  values  of 
fv  up  to  twice  the  values  given  in  Figs.  5-8. 
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A  number  of  calculations  have  been  performed 
with  the  small-disturbance  code  for  transonic 
flows  with  sinusoidal  and  sharp-edged  gusts,  whose 

distributions  of  vorticity  are  small  in  amplitude, 

1 2 

but  will ch  are  not  necessarily  irrotational . 

These  correspond  to  examples  that  have  previously 
appeared  in  the  literature  for  Incompressible  or 
subsonic  flows. 


Sinusoidal  Transverse  Gusts 

The  first  case  to  be  considered  is  the  tran¬ 
sonic  counterpart  of  the  well-known  sinusoidal 
traveling-wave  gust  with  a  transverse  velocity 

vr  *  A1  sin(wt  -  uix)  0*1) 

where  x/c  -  2x/c  -  1  is  measured  from  the  mid¬ 
chord  of  the  airfoil,  and  wherr.  o  is  the  chord. 
This  case  is  particularly  straightforward  within 


the  framework  of  the  transonic  small-disturbance 
approximation,  3lnce  the  gu3t  velocity  does  not 
appear  explicitly  in  the  governing  equation.  The 
effect  of  the  gust  is  felt  through  the  boundary 
condition  on  the  airfoil  as  an  effective  time-  and 
space-dependent  angle  of  attack,  Eq.  (12). 

A  sketch  of  this  problem  and  typical  results 
for  this  type  of  gust  are  shown  in  Fig.  9;  the 
figure  shows  computations  for  both  the  nonlinear 
and  linear  mode  of  Eq.  (11)  at  various  "phases" 
during  the  passage  of  the  gust.  This  figure  may 
be  compared  with  Fig.  5,  for  which  the  airfoil, 
Mach  number,  and  mean  angle  of  attack  are  the 
same.  Here  the  nondimensional  frequency  param¬ 
eter,  K  -  wc/U^,  is  inversely  proportional  to  the 
wave  length  A  of  the  gust;  X  -  2*11^/10  *  4mc  in 
this  particular  case.  The  amplitude  of  the  gust 
corresponds  to  an  effective  fluctuating  angle  of 
attack  of  ±1 0 . 


Fig.  9  Transonic  airfoil-gust  interaction  on  the  NACA  0012  airfoil  at  M^  ■  0.8,  a  -  0;  sinusoidal  gust 

with  A,  -  0.017,  uic/UB  -  0.5. 
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(b)  Lower  surface. 

Fig.  11  Harmonic  components  of  the  pressure  dis¬ 
tributions  in  supercritical  and  conventional  lir- 
foils  at  M  *  0.76;  sinusoidal  gust  with 
A,  >  0.017,  uc/’Jm  -  0.5. 

oscillating  airfoil  problems,  if  the  shock-wave 
motion  is  negligible  and  the  unsteady  fluctuations 
are  essentially  linear  perturbations  about  a  non¬ 
linear  steady  or  mean  ‘‘low.  In  this  section  we 
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discuss  the  indictal  response  to  a  sharp-edged 
gust.  Fig.  12,  in  subsonic  and  transonic  flow. 


Fig.  12  Schematic  of  indicial  response  for 
plunging  motion  (above)  and  for  a  sharp-edged  gust 
(below) . 


Figure  12  illustrates  the  similarities  and 
differences  in  the  response  of  a  thin  airfoil  to 
two  types  of  indicial  inputs;  namely, 

{0  for  t  <  0 

( lb) 

A1  for  t  >  0 

fO  for  Uat  <  X 

gust;  vc  *  J  (16) 

l  A,  for  U^t  >  X 

According  to  linear  theory,20  the  sharp-edged  gust 

response  rises  monotonically  froiji  zero  at  t  -  0 

to  the  value  C.  =  2itA,/(1  -  M2)  a3  t  ■»  ", 
LI" 

whereas  the  plunge  response  starts  at  -  nA1 
for  M  "0  and  at  9A./M  for  M  >  0.  In 

CO  I  00  30 

either  case,  the  approach  to  the  steady-state 
value  is  progressively  slower  as  increases. 

Figure  13  shows  the  nonlinear  gust  response 
of  the  NACA  69A006  airfoil  as  a  function  of  Mach 
number,  with  a  logarithmic  abscissa.  The  critical 


Mach  number  for  this  airfoil  is  approximately  0.83 
for  a  -  0,  and  the  flow  is  only  weakly  transonic 
at  -  0.80  and  a  -  1°.  Therefore,  the  lift 
response  shown  in  Fig.  13  is  essentially  linear 
for  M  <  0.80. 

<B  — 

However,  the  behavior  becomes  distinctly 
nonlinear  at  higher  Mach  numbers.  The  overshoot 
in  lift  at  M^  -  0.9  is  due  to  a  complex,  non¬ 
linear  motion  of  the  shock  wave  on  the  lower  sur¬ 
face,  as  shown  in  Fig.  19.  The  upper-surface 
shock  wave  moves  monotonically  rearward  from  its 
initial  position  of  x  «  0.85  and  stabilizes  at 
the  trailing  edge  at  U^t/c  »  25.  However,  the 
shock  wave  on  the  lower  surface  moves  forward 
initially,  reaches  a  minimum  value  of  Xg  at 
U^t/c  -  30,  and  then  undergoes  a  small  oscillation 
in  position  between  U^t/c  •  90  and  60,  before 
reaching  its  final  rearward  location.  This 
nonlinear  behavior  indicates  that  it  would  be  a 
mistake  to  use  the  Indicial  gust  response  at 
M^  -  0.9  in  a  superposition  integral  for 
arbitrary  wave  forms. 


.01  1  1  10  100  1000  10  000 


U»«/c 

Fig.  13  Indicial  gust  response  on  the 
NACA  69A0012  airfoil . 
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disturbances  superimposed  on  a  uniform  free 
stream,  under  the  assumption  that  the  prescribed 
disturbance  remains  undistorted  as  it  convects 
past  a  thin  airfoil.  A  number  of  unsteady  alr- 
foil-gust  interactions  have  been  calculated  within 
this  framework,  and  sample  results  illustrate  the 
essential  effects  of  periodic  gusts  on  the  tran¬ 
sonic  flow  around  conventional  and  supercritical 
airfoils.  Also,  several  examples  of  sharp-edged 
gusts  have  been  studied.  As  in  the  concentrated- 
vortex  interactions,  the  strength  and  unsteady 
motion  of  the  shock  wave  were  found  to  play  major 
roles  in  the  flow-field  development  and,  conse¬ 
quently,  in  the  airloads  on  the  airfoil. 
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:  16  Abstract 


Unsteady  interactions  ot  strong  concentrated  vortices,  distributed 
gusts,  and  sharp-edged  gusts  with  stationary  airfoils  have  been  analyzed  in 
two-dimensional  transonic  flow.  A  simple  and  efficient  method  for  introduc¬ 
ing  such  vortical  disturbances  has  been  implemented  in  numerical  codes  that 
range  from  inviscid  transonic  small  disturbance  to  thin-layer  Navi.er  Stokes. 
The  numerical  results  demonstrate  the  large  distortions  in  the  overall  flow 


io  and  in  the  stir  t  ace  air  !  e:d'.  I  li.it  are  produced  bv  various  vortical 
interact i eus.  The  results  oi  tile  d  i  1 fereot  codes  ale  in  excellent  qua  I  i  ta- 
t i ve  agreement ,  but,  as  might  expected,  the  transonic  sma 1 1-d i sturhanee 
calculations  are  deficient  in  tin  important  region  near  the  leading  edge. 
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